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Biomass Partitioning Affects the Growth of Pinus 
Species from Different Elevations 


ZHANG Shi-Bao 


( Key Laboratory of Economic Plants and Biotechnology, Kunming Institute of Botany , 
Chinese Academy of Sciences, Kunming 650201, China) 


Abstract; The conifer forests in southwestern China are mainly dominated by three vicariant species within Pinus : 
P. kesiya var. langbianensis, P. yunnanensis, and P. densata. Their sites range from lower to higher elevations, re- 
spectively, and each species shows differences in rates of development, especially with regard to height. To identify 
the physiological and morphological causes of this inherent variation, photosynthesis, biomass partitioning, growth 
rates and leaf traits were investigated of plants cultivated under the same environmental conditions. Trees of the spe- 
cies native to the lower elevation were taller, and had higher values for dry weight, relative growth rate (RGR) , leaf 
mass fraction ( LMF) , stem mass fraction (SMF), and specific leaf area per unit mass (SLA), relative to those 
from the higher elevations. However, their leaf N and C contents per unit area, and their root mass fraction (RMF) , 
were smaller than those of high-elevation trees. Photosynthetic capacity in species from high elevations was not signif- 
icantly reduced from the level calculated for trees from lower elevations. Both RGR and tree height were positively 
correlated with LMF and negatively with RMF, but no significant positive correlations were found with maximum pho- 
) and mass-basis (A 


tosynthetic rate determined on both an area-basis (A ). These findings suggest that the 


max mass 


patterns of biomass partitioning and long-term morphological traits are better predictors of performance among trees of 
different Pinus species growing along an elevational gradient. 
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Wide variations in growth rates among trees 
across species are closely correlated with investments 
of carbon resources and ecological distributions 
(Wright and Westoby, 2000; Poorter and Garnier, 
2007; King et al., 2013). Fast-growing species are 
generally found in relatively more fertile habitats , 
whereas species that occupy infertile environments 
tend to have low growth rates ( Biere, 1996; Poorter 
and Garnier, 2007). Those rates can affect seedling 
survival, reproduction, productivity, competition, 
and forest structures (Lambers and Poorter, 1992; 
Biere, 1996; Coomes and Allen, 2007). Therefore , 
examining the factors that influence tree growth is 
essential for modeling forest productivity and its 
functioning under climatic change (Coomes and Al- 
len, 2007; King et al., 2013). 

Relative growth rate (RGR) is a complex trait 
that is determined by differences in physiology , mor- 
phology, and biomass partitioning (Shipley, 2006; 
Poorter and Garnier, 2007; Poorter et al., 2012). 
This trait encompasses three variable components: 
net assimilation rate (NAR), specific leaf area 
(SLA), and leaf mass ratio (LMR) ( Poorter and 
Garnier, 2007; Poorter et al., 2012; Tomlinson et 
al., 2012). For example, species from humid envi- 
ronments partition more biomass to the roots and less 
to the stems than those from semiarid environments 
in Australia, Africa, and South America (Tomlinson 
et al., 2012). Poorter and Remkes (1990) have 
found, in a controlled experiment, that the RGRs of 
69 plant species are most strongly correlated with 
SLA and LMR, but are also positively associated 
with biomass allocation to the leaves and NAR. How- 
ever, based on the results of a meta-analysis with 
614 species from 83 different experiments in Eu- 
rope, America, and Australia, Shipley (2006) has 
suggested that LMR is never strongly related to RGR 
but, instead, NAR is the best general predictor of 
variations in interspecific RGRs. Consequently, the 
relative contribution of SLA, NAR, and LMR to 
RGR varies among species or when plants are grown 


in contrasting environments. 


Trees in the Pinus genus cover a large geographi- 
cal area and a wide variety of habitats in southwes- 
tern China, where they are of great economic and ec- 
ological importance. As the dominant component of 
conifer forests in that region, the three vicariant spe- 
cies, from low to high elevation, are P. kesiya Royle 
ex Gordon var. langbianensis Gaussen , P. yunnanensis 
Franch., and P. densata Mast. (Wu, 1990). These 
species exhibit different rates of development in their 
habitats, especially in height increment, with P. 
kesiya being the fastest grower ( Dai et al., 2006). 
However, the physiological and morphological causes 
underlying those differences are not completely known. 

Plant heights and biomass production in trees 
generally decline with increasing elevation, as shown 
in both common garden studies and natural forests 
(Oleksyn et al., 1998; Li et al., 2003; Coomes and 
Allen , 2007). However, the extent of the correlation 
between elevation and annual stem elongation can dif- 
fer among species ( Angert, 2006). For example, 
highland plants might display slower growth rates be- 
cause of morphological or physiological reasons ( At- 
kin et al., 1996a; Oleksyn et al., 1998; Angert, 
2006). Hoch et al. (2002) have suggested that a low 
temperature-driven sink explains the lag in growth by 
P. cembra at higher elevations at the tree line in the 
Swiss Alps. Although some studies demonstrated that 
leaf thickness, leaf N content, and photosynthetic 
rates are greater in highland plants than in lowland 
plants on the island of Hawaii or the Northern Ameri- 
ca continent (Cordell et al., 1999; Hultine and Mar- 
shall, 2000) , conflicting results from other investiga- 
tions shown that a general trend in photosynthetic ca- 
pacity does not occur across elevations in the tropical 
high Andes and southwest China ( Cabrera et al., 
1998; Zhang et al., 2007; 2011). For example, 
clones of spruce from the tree line have a 4.3-fold 
lower growth rate and they contain 60% less chloro- 
phyll per unit mass than trees from valleys in Europe 
(Polle et al., 1999). Westbeek et al. (1999) have 
reported, from a common garden, that RGR is nega- 


tively correlated with leaf N content per unit area, as 
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well as with chlorophyll and Rubisco contents, for 
Poa species in both alpine and lowland regions. De- 
spite greater photosynthetic rates in high-elevation 
populations, seedling heights and dry masses meas- 
ured in common garden studies decline with the ele- 
vation at which their seed originated. Proportional dry 
mass partitioning in the roots nearly doubles with in- 
creasing elevation of origin in controlled experiments 
(Oleksyn et al., 1998). Consequently, researchers 
still do not have a clear understanding of the causes 
for low relative growth rates by highland plants. 
Here, the growth rates, leaf N contents, spe- 
cific leaf areas, biomass partitioning, and photosyn- 
thetic rates of P. kesiya, P. yunnanensis and P. den- 
sata were monitored in a common garden at Kunming 
Botanical Garden, China. The objective was to iden- 
tify the factors affecting the development of pines 
that originate from different elevations. We hypothe- 
sized that there is an inherent difference in the rela- 
tive growth rate of these three species, and that this 
difference is correlated with leaf physiology and bio- 


mass partitioning. 


1 Materials and methods 
1.1 Plant species and sites 

Three species in Pinus—P. kesiya, P. yunnanen- 
sis, and P. densata—were selected. All are evergreen 
species that dominate conifer forests, from low to 
high elevation, respectively, in southwestern China 
(Table 1). Seeds of P. kesiya, P. yunnanensis, and 
P. densata were collected from Puwen (elev. 1 377 
m), Jianchuan (2 196m), and Deqing (3 447 m), 
respectively, from October to December in 2006. 
They were sown at Kunming Botanical Garden (elev. 
1990 m; E 102. 74°, N 25. 15°) on 22 March 2007 
and had germinated by 9 April. 70 seedlings per 


species were transplanted on 7 May into plastic pots 
filled with a peat, forest soil, and humus mixture 
(1:3:1; vi viv). The initial dry mass of sample 
seedlings from each species was measured 1 d before 
transplanting occurred, and the final harvest was 
made on 12 October 2007. After transplanting, the 
seedlings were grown under full sunlight, and were 
fertilized monthly with a liquid nutrient solution and 
watered every 2 to 3 d. From May to October, the mean 
monthly temperature ranged from 13.2 °C to 19.9 °C 
(mean + SD, 17.6+2.4 °C) while mean monthly pre- 
cipitation was 19.6 mm to 204.0 mm (mean + SD, 
115.7+74.6 mm). 
1.2 Physiological measurements 

Gas exchange in response to light and CO, con- 
centration was recorded on 7 to 9 October 2007 from 
fully expanded leaves, using a LI-6400 portable 
photosynthesis system (LI-COR, Lincoln, NE, USA). 
Data were recorded from 10 plants per species. Be- 
fore the measurements began, each sample leaf was 
exposed to actinic light of 1 200 umol m~s' (10% 
blue, 90% red) for 15 min to induce maximum 
stomatal opening. Curves for the photosynthetic light 
response (A-PPFD) were made using an automated 
protocol built into the LI-6400. The program was 
configured to advance to the next step if the sum of 
three coefficients of variation (CO,, water vapor, 
and flow rate) was less than 0.3%. The minimum 
waiting time was 3 min. Each leaf was equilibrated to 
initial conditions by waiting at least 15 min before 
executing the automated protocol. A-PPFD curves 
were generated for 6 leaves per species, at light in- 
tensities of 2 000, 1 600, 1 200, 1000, 800, 600, 
400, 300, 200, 100, 50, and 0 pmol m~s'. Other 
test conditions included a controlled level of CO, 


(380 pmol mol!) , a flow rate of 500 pmol s™', a 


Table 1 Origins of the three Pinus species used in this study 


Site for seed collection 








Species Code Elevation range/m 
Longitude Latitude Elevation/m 
P. kesiya Pk 600-1600 E100°55' N22°22’ 1377 
P. yunnanensis Py 1000-3200 E99°54' N26°32’ 2196 
P. densata Pd 2600-3600 E98°52’ N28°27’ 3447 
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leaf temperature of 25 °C , and a vapor pressure defi- 
cit ( VPD) of 0. 6 to 1.0 kPa. 

Photosynthetic CO, response curves (A-C,) and 
A-PPFD curves were determined for the same leaves. 
After measurements were completed for the A-PPFD 
curve, each leaf was exposed for 15 min to a light 


1 


intensity of 1 200 pmol ms! and a CO, concentra- 
tion of 380 umol mol™'. Other conditions were a leaf 
temperature of 25 °C and VPD of 1.0 to 1.5 kPa. The 
A-C, curve measurement was started at ambient CO, 
concentration, which decreased gradually to 0 wmol 
mol ', returned to 380 umol mol’, and then in- 
creased to a higher concentration to ensure that the 
stomata remained open throughout the recording pe- 
riod. Photosynthetic rates were measured at different 
CO, concentrations, using the automated protocol 
built into the LI-6400. 

After completing the gas exchange measure- 
ments, we recalculated the photosynthetic rate based 


on actual leaf area values, which were estimated by 


the method of Johnson (1984) : 


LA = 2L(1 + m/n) /(nv/TL) (1) 


where LA is leaf area (cm~), L is needle length 
(cm), n is the number of needles per fascicle (n= 
3 for P. kesiya and P. yunnanensis; n=5 for P. den- 
sata) , and needle volume (v) is estimated via water 
displacement (Johnson, 1984). 

Photosynthetic light response curves were fitted 
with non-rectangular hyperbola. The maximum pho- 


tosynthetic rate (Aa) and respiration rate (R, ) 


were determined via Photosyn Assistant v1. 1 ( Dund- 
ee Scientific, Dundee, Scotland, UK) according to 
the method of Prioul and Chartier (1977). Using A- 


C, curves, we calculated the maximum carboxylation 


rate by Rubisco ( V 


cemax 


) and light-saturated electron 


transport (Jax) with Photosyn Assistant, based on 


the photosynthetic model of von Caemmerer and Far- 
quhar (1981). 

Mesophyll diffusion conductance (g,,) from the 
internal air space to the chloroplasts was estimated 


according to the method of Harley et al. (1992) ; 


A 
&m T * (2) 
r Kaa +8(A +Rı)] 
i Tiis = 4(A + Ra) 


where the value for R, is found from the A-PPFD 








curve, and I” is the hypothetical CO, compensation 
point in the absence of R} (42.75 pmol mol” at 25 °C) 
( Bernacchi et al., 2001). Values for g,, were calcu- 
lated from our measurements of photosynthesis at in- 
ternal CO, concentrations of 100 to 300 pmol mol”, 
with the average value of g„ being determined for 
each leaf. 

Chlorophyll was extracted per the technique of 
Moran and Porath (1980), and its concentration 
was analyzed on a UV-2550 spectrophotometer ( Shi- 
madzu, Japan) before being calculated according, to 
the equations of Inskeep and Bloom (1985). 

1.3 Growth rate and biomass partitioning 

After the photosynthetic data were obtained, the 
plants were harvested and divided into root, stem, 
and leaf portions. Ten plants were measured for each 
species. The dry mass of each portion was recorded 
after the tissues were dried for 48 h at 80 °C. After- 
ward, the leaf mass fraction (LMF) , stem mass frac- 
tion (SMF), and root mass fraction (RMF) were de- 
termined in proportion to total dry mass per plant. 
The values for relative growth rate were calculated by 
the following formula (Poorter and Garnier, 2007). 

InM, - InM, 


t =t 


RGR (3) 


where M, and M, are the biomass at time ¢, and t,, 
respectively; and ż, -t is the time span between 
measuring events. Here, that span was 185 d. 

Leaf area per plant was determined according to 
the method of Johnson (1984). From this, SLA was 
calculated as the specific leaf area per unit mass 
(m? kg”). Leaf N content was assessed with a Leco 
FP-428 nitrogen analyzer ( Leco Corporation, St. Jo- 
seph, MI, USA). 

1.4 Statistical analysis 

Statistical analysis was performed with SPSS 13.0 
(SPSS Inc., Chicago, IL, USA). We used one-way 
ANOVA and LSD multiple comparison tests to esti- 
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mate the differences in RGR, biomass partitioning , 
physiological parameters, and leaf N content. Rela- 
tionships between RGR and leaf traits were assessed 
through Pearson’ s regression analysis. A principal 
component analysis (PCA) was performed to char- 


acterize the associations among 17 leaf traits. 


2 Results 
2.1 Interspecific variations in growth, biomass 
partitioning , and physiology 

Tree height (H) and dry mass per plant varied 
significantly across species (Fig. 1). For example, 
plants of P. kesiya were taller than those of P. yun- 
nanensis and P. densata ( P<0. 001) whereas the mean 
basal diameter from P. yunnanensis (1. 179+0. 140 cm) 
was larger than from P. densata (0. 747+0. 073 cm). 
However, the dry mass per plant and the RGRs of 
P. kesiya and P. yunnanensis were greater than those 
of P. densata (Fig. 1). Consequently, species from 


the lower elevations grew more rapidly than the one 
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Fig. 1 Differences in plant height, dry mass weight (DW) , relative 
growth rate (RGR), leaf mass fraction (LMF) , stem mass fraction 
(SMF), and root mass fraction (RMF) for Pinus kesiya (Pk), P. 
yunnanensis (Py), and P. densata (Pd) grown under the same envi- 
ronmental conditions. Different letters above bars for each component 
indicate statistically significant differences in mean values (P <0. 05) , 


as determined by LSD multiple comparison tests 


from the highest elevation. 

Biomass partitioning to the organs also differed 
significantly among species (Fig. 1). Whereas the frac- 
tions of dry matter produced in the leaf and stem were 
higher in P. kesiya and P. yunnanensis than in P. densa- 
ta, the root mass fraction showed an opposite trend. 

Significant differences among species were found 
for the area-based respiration rate (R,) and stomatal 
conductance (g,), but not for the area-based photo- 
synthetic rate (A,,,, ), mass-based photosynthetic 
rate (A 


maximum carboxylation rate ( V 


emax 


), mass-based respiration rate (Ramas), 


) , light-saturated 
electron transport rate (Jaa), and mesophyll con- 


ductance (g,,) (Fig. 2). 


max 
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Fig.2 Comparisons of area-based photosynthetic rate (Anax) , mass- 
based photosynthetic rate ( Amas ), maximum carboxylation rate 


CV ciatie) area-based 


, light-saturated electron transport rate (Jia. ); 
respiration rate ( R4) , mass-based respiration rate (Ry jac.) stomatal 
conductance (g,) , and mesophyll conductance (g,, ) among Pinus 
kesiya (Pk), P. yunnanensis (Py), and P. densata (Pd). Different 
letters above bars for each component indicate statistically significant 


differences in mean values (PSO. 05) , as determined by LSD multi- 


ple comparison tests 
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2.2 Determinants of growth rate 

Tree height was correlated significantly and posi- 
tively with LMF, SMF, SLA and RGR, but nega- 
tively with RMF (Fig. 4). No significant correlation 
existed between H and Ana: 

Dry weight per plant was correlated positively 
with LMF but negatively with RMF. Neither R, nor 
Aa Was correlated with dry weight (Fig. 5). RGR 
was significantly and positively correlated with LMF 
but negatively with RMF. No significant correlations 
were found between RGR and Ana or Anass (Fig. 6). 
Finally, A,,,, was positively correlated with Vmax 
Jax Zax Zn, and leaf N content per unit area, but 


max 


negatively with SLA (data not shown). 
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Fig.5 Pearson correlations of dry mass weight (DW) with leaf mass 
fraction (LMF) , root mass fraction ( RMF) , maximum photosynthetic 
rate (A 


max)» and respiration rate (R,) for 3 Pinus species 


Principal component analysis showed that H, 
RGR, LMF, RSF, RMF, and SLA loaded mainly 
on the first PCA axis, explaining 36. 6% of the total 
variation; Á maxs Jimax > 


V max ’ Es ’ Em ’ and leaf N con- 


tent per unit area (Na) loaded on the second axis, 


explaining 26. 3% of the total (Fig. 7). The first axis 
of the PCA was mainly associated with plant growth 
and biomass partitioning while the second axis was 


associated with photosynthetic carbon assimilation. 
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Fig.6 Pearson correlations of relative growth rate ( RGR) with leaf 
mass fraction (LMF) , root mass fraction ( RMF) , mass-based maxi- 
mum photosynthetic rate (Amass ) , and area-based maximum photosyn- 
thetic rate (A 


) for 3 Pinus species 
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Fig.7 Principle component analysis of 17 leaf traits from 3 Pinus 


species. A,,,.,, gross photosynthetic rate; A area-based maximum 


max» 
photosynthetic rate; DW, dry mass weight per plant; g,,, mesophyll 
conductance; g,, stomatal conductance; H, tree height; J,,,,, light- 
saturated electron transport rate; LMF, leaf mass fraction; N 


area » Nl- 


trogen content per unit area; R}, dark-respiration rate; RGR, relative 
growth rate; RMF, root mass fraction; R,, photorespiration rate; 
RSF, root-shoot ratio; SLA, specific leaf area per unit mass; SMF, 


stem mass fraction; V, 


emax ? 


maximum carboxylation rate 


3 Discussion 
3.1 Growth rate in relation to elevation 
Seedlings of Pinus species originating from the 
highest elevation were shorter and had smaller values 
for RGR than congeneric species from the lowest ele- 
vation. This indicated that the former grew more 
slowly. These findings are in line with previous re- 
ports from common gardens and natural forests that 
tree height and biomass production are reduced as 
elevation increases ( Oleksyn et al., 1998; Li et al., 
2004; Coomes and Allen, 2007). Li et al. (2003) 


have found that the mean annual biomass increment 


per tree in the eastern Himalayas is less as elevation 
increases, i.e., reduced by 476 g per 100 m between 
1 680 and 1 810m and by 103 g between 1 810 and 
1 940 m. However, Angert (2006) has shown in 
common garden experiments that trees within a given 
species accumulate their greatest aboveground bio- 
mass when grown under a temperature regime that is 
characteristic of the center of its natural elevation 
range. Thus, the effect of seed origin or elevation on 
RGR is considered species-specific ( Angert, 2006). 
3.2 Growth rate in relation to biomass partitioning 

The decline in growth rate with increasing eleva- 
tion can be caused by both environmental and genetic 
factors (Korner, 2003; Li et al., 2003; King et al., 
2013). For example, in natural forests, a shorter 
growing season that results from low temperatures at a 
higher elevation can reduce a tree’s accumulation of 
C and its growth rate ( Oleksyn et al., 1998; Körner, 
2003). Because seeds of the trees used in our study 
were collected at different elevations but were then 
exposed to the same experimental conditions, any 
interspecific variations in tree growth rates would 
have reflected inherent genetic differences. 

Present study demonstrated that tree height, dry 
weight, and relative growth rate were correlated posi- 
tively with the leaf mass fraction, but negatively with 
the root mass fraction. Previous studies under con- 
trolled conditions have also revealed that RGR is pos- 
itively associated with biomass partitioning in the 
leaves (Poorter and Remkes, 1990; Wright and Westo- 
by, 2000). By contrast, Shipley (2006) has noted 
no significant correlation between RGR and LMF in 
83 different experiments. The essential resources 
needed by plants are acquired by different organs; 
more resources allocated to the leaves can increase 
light capture, while more resources transported to 
the roots can improve the uptake of water and miner- 
al elements (Taub, 2004). Because plants are com- 
posed largely of materials derived via photosynthesis , 
an increased partitioning to non-photosynthetic organs 
would reduce their growth rate in general ( Mooney , 


1972). Consequently, the functional balance between 
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leaves and roots is an important determinant of inter- 
specific differences in RGR (Wright and Westoby, 
2000; Osone et al., 2008). 

Biomass partitioning in different organs may re- 
flect a survival strategy that is adjusted to meet cur- 
rent environmental conditions ( Poorter et al., 2012). 
The process of allocating C to either the roots or 
shoots is more sensitive to temperature, shade, and 
water deficit ( Poorter et al., 2012; Tomlinson et al., 
2012). Generally, plants in cold environments tend 
to accumulate more dry matter in the roots, boosting 
their production of those tissues ( Oleksyn et al., 
1998; Li et al., 2004), while also reducing their 
stem and leaf fractions to increase RMF values 
( Poorter et al., 2012). Plants that decrease the pro- 
portion of leaves and partition more biomass to their 
roots at higher elevations or on less productive sites 
will benefit from faster recovery of their growth rates 
the following year, and will be better able to adapt to 
unfavorable climatic conditions ( Oleksyn et al., 
1998; Korner, 2003). However, improved survivabil- 
ity comes at the expense of slower growth because 
plants must respond to environmental gradients by ad- 
justing their pattern of partitioning biomass to maxi- 
mize growth performance ( Wang et al., 2008). 

3.3 Correlation of growth rate with photosyn- 
thetic rate 

Maximum photosynthetic rate was not correlated 
with tree height, plant dry weight, or RGR in our 
three Pinus species. These results contradict previous 
findings from controlled experiments that RGR is cor- 
related with the photosynthetic rate ( Lambers and 
Poorter, 1992; Atkin et al., 1996b; Shipley, 2006). 
Villar et al. (2005) have also suggested that RGR is 
correlated with NAR during a short growing period, 
but with morphological traits over the long term. How- 
ever, maximum photosynthetic rates for the three spe- 
cies did not vary with elevation, but were positively 
correlated with CO, diffusive conductance and leaf N 
content per unit area while negatively correlated with 
SLA. Earlier researchers have implied that no gener- 


al trend in photosynthetic rate and leaf traits exists 


across elevations in natural habitats (Cabrera et al. , 
1998, Zhang et al., 2007; 2011). Villar et al. (1998) 
have suggested that the photosynthetic rate in Ae- 
gilops species is negatively correlated with SLA in a 
common garden, thereby accounting for the absence 
of a correlation between SLA and RGR. 

In conclusion, the pines grown from seed gath- 
ered at higher elevations had slower growth rates, 
but their photosynthetic capacity was not significantly 
inferior to that of trees representing lower elevations. 
We believe that, at higher elevations, plants tend to 
allocate more biomass to the roots while reducing 
their leaf production, which ultimately leads to slo- 
wer growth. The pattern of biomass partitioning, 
rather than photosynthetic capacity, determines any 
interspecific variations in growth rates among Pinus 
species originating from different elevations. Thus, 
long-term morphological traits are better predictors of 


growth rates for trees along an elevational gradient. 
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